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I.  BASIC  DESCRIPTION 


The  computer  code  (Version  2.0  of  SMLJi.O)  described  in  this  technical 
report  was  developed  as  part  of  research  projects  under  the  direction  of 
Dr.  A.  Chames  (The  University  of  Texas  at  Austin).  Office  of  Naval 
Research  Contract  N00014-80-C-0242,  and  Dr.  W.W.  Cooper  (The  University 
of  Texas  at  Austin),  Office  of  Naval  Research  Contract  N00014-76-C-0932. 
Dr.  Authella  Bessent  (The  University  of  Texas  at  Austin),  Dr.  Wailand 
Bessent  (The  University  of  Texas  at  Austin),  and  Dr.  J.  Kennington 
(Southern  Methodist  University)  joined  the  group  during  1979  and 
designed  and  tested  SMU1.0  during  the  last  six  months  of  that  year. 

Given  a  set  of  decision  making  units  (DMU's)  with  common  inputs 
and  outputs,  SMU1.0  will  compute  the  efficiency  of  each  of  the  DMU's. 

The  underlying  theory  is  described  in  Chames,  Cooper,  and  Rhodes  [1]. 

S MU 1.0  is  written  in  FORTRAN  and  was  tested  on  SMU's  CDC  Cyber  73 
using  the  MNF  compiler.  The  code  employs  a  modified  primal  simplex  code 
which  was  originally  developed  by  Iqbal  All  (The  University  of  Texas  at 
Austin)  as  a  part  of  his  dissertation  work.  The  input  format  and  output 
reports  were  designed  by  Authella  and  Wailand  Bessent. 


II.  MODEL  AND  SOFTWARE  DESIGN 


To  determine  the  efficiency  of  any  DMU  requires  t;he  solution  of  a 

particular  linear  program.  In  this  section,  we  present  the  mathematical 

model,  we  give  the  initial  feasible  basis  used  to  initiate  the  primal 

simplex  algorithm,  and  we  briefly  discuss  the  design  of  the  system. 

A.  Notation 

m  "  number  of  DMU's 

n^  ■»  number  of  outputs 

n2  “  number  of  inputs 

A  »  an  n^  x  m  matrix  of  outputs 

C  *  an  n2  x  m  matrix  of  inputs 

A(j)  ■  the  jth  column  of  A 

C(j)  =  the  jtl1  column  of  C 

a^j  ■  the  (i,  j)th  element  of  A 

cij  *  the  (i,  j)1"*1  element  Of  C 

e.  "a  column  vector  having  i^  component  equal 
+1  and  all  other  components  zero. 


Assumptions 


a^  i  0  for  all  i,  j. 
cij  i  0  for  all  i,  j. 
C.  Model 


Consider  the  following  linear  program: 


X 


where  X  is  an  m  vector,  sx  is  an  nj  vector  of  slacks,  and  s2  is  an  n2 
vector  of  slacks.  Let  (a*,  X*,  s*,  s*)  denote  an  optimum  for  P(j). 


If 


z  >1,  then  DMU  j  is  said  to  be  efficient; 
* 


z  >1,  then  DMU  j  is  said  to  be  inefficient. 

D.  Initial  Feasible  Basis 

The  following  matrix  is  a  feasible  basis  for  P(j)  with  *  ■  X^  ■  1 
and  all  other  variables  equal  zero: 


r  l 
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E.  Software  Design 

SMU1.0  solves  P(l)  through  P(m)  via  the  revised  primal  simplex  method 
using  the  starting  bases  (1)  .for  eaqfo  problem,  P(j).  The  basis  inverse  is 
maintained  in  product  form  and  only  the  non-zeroes  are  stored.  The  reinver¬ 
sion  routine  is  based  on  the  work  of  Hellerman  and  Rarick  [4].  This  routine 
incorporates  a  technique  known  as  "splitting  the  bump"  (see  Kalan  [5]  and 
Orchard-Hays  [6])  and  uses  the  "spike  swapping  theory"  of  Helgason  and 
Kennington  [2].  A  good  description  of  the  computational  procedures  may  be 
found  in  Helgason  [3].  The  arrays  and  working  files  require  approximately 
(3m  +  13)  (n^  +  n2)  +  4m  words  of  core  storage. 
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III.  INPUT  FORMAT 


Card  Group  Source  File  Composition 


Description  and  Format 


1  Tape  1  1  card 

2  Tape  1  1  card 


3 


Tape  1 


card  group 


4 


Tape  1 


1  card 


"Problem  Title"  (13A6,  A2) 

"Number  of  DMU's",  "Number  of 
Outputs",  "Number  of  Inputs" 
(315) 

"Output  1  Label",  "Output  2 
Label",  ...,  "Output  Last 
Label",  "Input  1  Label", 
"Input  2  Label",  . . . ,  "Input 
Last  Label"  (8A10) 

"Blank  Card" 


5 


Tape  1 


card  group  Outputs  for  the  DMU  followed 

for  each  DMU  by  Inputs  for  the  DMU  (I0F8.0) 


Card  groups  6  and  7  are  optional.  If  included,  each  DMJ's  name  will 
appear  on  its  output  report. 


6 

7 


Tape  I 


Tape  1 


1  card  "Blank  Card"  (Use  only  if 

card  Group  7  follows) 


card  group  "DMU1  name"  (up  to  80  characters) 

for  each  DMU  "DMJ2  name" 

(in  same  order  ' 

as  DMU  data)  ' 

t 

"DMUm  name" 


Card  group  8  is  optional.  A  blank  card  defaults  to  an  analysis  of  all  DMU's. 


8  Input  1  card  "First  DMU  to  be  evaluated  in 

this  run",  "Last  DMU  to  be 
evaluated  in  this  run"  (free 
format,  i.e.,  the  two  nuinbers 
may  be  in  any  columns  as  long 
as  they  are  in  the  correct 
sequence  and  separated  by  either 
a  blank  or  a  comma) 
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IV.  SET  CODE  DIMENSIONS 

To  sec  Che  dimension  for  a  problem  one  should  use  a  CexC  edicor  which 
has  Che  capablliCy  of  replacing  all  occurrences  of  a  specified  characcer 
sCring  by  a  number. 

A.  The  Main  Working  SCorage 

Replace  "ZMAT"  wich  " - ".  (should  be  aC  lease  (m  +  nj  +  n2  +  3) 

(nj  +  n2  +  2)  +  (nj  +  n2  +  2) 

B.  For  The  Number  of  InefficienC  DMU's 

Replace  "ZINEFF"  with  (ra  is  sufficienCly  large) 

C.  For  Column  Lengch  Arrays 

Replace  "ZCOL"  with  "- - (Should  be  at  least  m  +  +  n2  +  3) 

D.  For  Row  Length  Arrays 

Replace  "ZROW"  with  " - ".  (Should  be  at  least  ^  +  n2  +  2) 

E.  For  Hie  Length  Of  The  ETA  File 

Replace  "ZETA"  with"-—-".  (Should  be  at  least  2m) 

WARNING:  For  each  of  the  above  replacements,  replace  the  words  with  a  number 
having  the  same  number  of  digits.  Use  leading  zeroes  if  necessary.  Suppose 
m  -  100.  Then  replace  ZETA  with  0200  and  replace  ZINEFF  with  000100. 
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